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ABSTRACT 


The  Free  Molecule  Micro-Resistojet  wsfs  designed  as  a 
micropropulsion  system  capable  /of  performing 
attitude  control  and  primary  / maneuvers  for 
nano  spacecraft  with^nassjjess  than  *10  kg.  The  details 
of  gas-surface  interactions  between  propellant 
molecules  and  surfaces  held  at  elevated  temperature 
are  critical  in  predicting  the  propulsion  system’s 
performance  and  efficiency.  The  aim  of  this  study  is 
to  asseST^arametrically  the  performance  of  a  typical 
thruster  geometry  using  a  general  Maxwell  scattering 
model  and  two  versions  of  the  Cercignani-Lampis- 
Lord  model.  The  models  are  incorporated  into  a 
Direct  Simulation  Monte  Carlo  numerical  code  and 
are  used  to  bound  the  predicted  performance 
characteristics  of  the  thruster.  The  total  specific 
impulse  varies  by  approximately  20%  over  range  of 
accommodation  coefficients  from  specular  to  diffuse 
surface  scattering.  However,  there  was  only  a 
maximum  difference  of  about  5%  between  the  models 
for  a  given  accommodation  coefficient.  Other  more 
microscopic  parameters,  such  as  axial  velocity 
distribution  functions,  appear  to  depend  more  on  the 
scattering  model  used. 
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INTRODUCTION 

VThe  Free  Molecule  Micro-Resistojet  (FMMR) 
has  been  designed  and  developed  as  a  primary 
and  attitude  control  thruster  for  nanospacecraft 
(m  <  10  kg).1*2  The  FMMR  is  constructed  using 
microelectromechanical  systems  (MEMS) 
fabrication  techniques  which  allows  it  to  be 
machined  on  a  micrometer  scale  and  to  be  easily 
integrated  with  other  MEMS  components  such  as 
embedded  control  systems,  valves,  and  pressure 
regulators.  Although  MEMS  fabrication  can  be 
very  beneficial  to  micropropulsion  systems,  a 
serious  drawback  of  traditional  MEMS 
micromachining  is  that  only  compatible 
materials  can  be  used  in  the  fabrication  process 
which  can  limit  device  performance.^  A 
schematic  of  the  FMMR  is  illustrated  in  Fi^y  1. 

The  propellant  gas  enters  into  the  FMMR 
stagnation  region  through  a  MEMS  valve  and 
filter  assembly  connected  by  the  propellant 
inlets.  The  device  operates  at  unusually  low 
pressures  (50  -  500  Pa)  which  gives  the 
propellant  gas  molecules  a  mean  free  path  on  the 
order  of  the  expansion  slot  width.  With  MEMS 
fabrication  techniques,  the  slot  size  can  be 
machined  with  a  width  as  small  as  -  1  pm; 
however,  the  nominal  FMMRpgeometrv  used  in 
this  study  has  a  slot  width  yi^TooTim.  Energy 
is  imparted  to  the  propellant  gas  through 
collisions  with  a  thin  film  heating  element 
suspended  above  a  silicon  pedestal.  The 
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pedestal  actjjto  minimize  direct  line-of-sight 
propellant  transport  through  an  expansion  slot 
after  entering  the  stagnation  region. 

The  design  requirement  is  to  arrange  a  surface 
held  at  the  required  temperature  to  be  the  last 
surface  contacted  by  a  propellant  molecule 
before  it  exits  the  device  through  an  expansion 
slot.  This  requirement  suggests  that  the  spacing 
between  the  thin  film  heater  and  the  expansion 
slot  be  on  the  order  of  the  molecular  mean  free 
path  in  the  device  to  reduce  intermolecular 
collisions  before  expansion. 

The  FMMR  relies  on  the  transfer  of  energy  into 
a  propellant  gas  through  molecular  collisions 
with  the  heated  surface.  Since  the  thin  film 
heater  is  the  only  heated  surface  in  the  device, 
subsequent  molecular  collisions  with  relatively 
cold  thruster  surfaces  is  expected  to  be  a  major 
loss  mechanism.  For  these  reasons,  the  details  of 
the  gas-surface  interaction  are  critical  in 
predicting  the  propulsion  system’s  performance 
and  efficiency.  The  aim  of  this  study  is  to  assess 
parametrically  the  performance  (thrust  and 
specific  impulse)  of  a  typical  FMMR  geometry 
using  several  gas-surface  interaction  models. 
The  models  are  incorporated  into  a  Direct 
Simulation  Monte  Carlo  (DSMC)3  code  to  bound 
the  predicted  performance  characteristics  of  the 
thruster. 


In  the  case  of  the  FMMR,  the  stagnation 
temperature  is  in  the  vicinity  near  the  expansion 
slot,  and  T*  is  the  temperature  of  the  molecules  at 
the  propellant  inlet  of  the  thruster.  The  free 
molecule  mass  flow  through  the  orifice  or  slot  is 


m  =  m 


(3) 


where  m  is  the  propellant  molecular  mass,  and 
T  is  the  average  thermal  speed  of  the  gas. 

A  propulsion  system’s  efficiency  is  generally 
evaluated  based  on  the  thruster’s  intrinsic 
specific  impulse  or  Isp.  The  Isp  of  a  propulsion 
system  is  the  ratio  of  the  thrust  generated  to  the 
propellant  mass  flow  required  to  produce  the 
thrust.  For  a  free  molecule  flow  through  a  thin 
orifice  or  slot  the  specific  impulse  is  given  by 


*kT0 
j  _  V  2m 


(4) 


where  g0  is  the  gravitational  constant.  As 
evident  in  the  above  formulation,  the  details  of 
the  energy  exchange  between  the  thin  film 
heaters  and  propellant  molecules  drive  the 
performance  of  the  FMMR. 


THEORY 

In  a  uniform  free  molecule  flow,  the  thrust  3,  or 
flux  of  normal  momentum,  produced  by  a  gas  at 
a  stagnation  pressure  of  p0  through  a  thin  orifice 
or  slot  of  area  A0  is 


where  n0  and  T0  are  the  number  density  and 
temperature  in  the  stagnation  region.  In  the  case 
of  the  FMMR,  the  stagnation  temperature  in  the 
vicinity  of  the  expansion  slot  is  governed  by  the 
physics  of  the  interaction  between  the  propellant 
gas  and  the  thin  film  heated  element.  For  the 
present  study,  the  thin  film  heated  element  is 
assumed  to  be  an  inert  surface  material  defined 
by  a  wall  temperature  Tw  and  a  single  energy 
accommodation  coefficient  a.  Therefore, 

(2) 


CALCULATIONS 

The  analytical,  free  molecule  results  from  the 
previous  section  are  expected  to  be  useful  in  the 
basic  design  of  the  FMMR;  however,  Direei 
Simulation  Monte  calculations 

are  required  to  predict  flowfield  properties  of  the 
actual  FMMR  under  transitional  rarefied  gas 
conditions.  Many  of  the  input  models  required, 
such  as  gas-surface  interaction,  are  topics  of 
current  research.  Parametric  variation  of  the 
model  features  can  be  used  to  assess 
performance  sensitivity  to  these  uncertainties 
independently. 

The  two-dimensional  computational  domain 
the  nominal  FMMR  geometry  is  shown  in  Fig/  2. 
The  geometry  is  symmetrical  about  the  slot 
centerline  IJ.  Boundary  AB  is  an  open  inlet; 
molecules  entering  are  selected  from  a  drifting 
equilibrium  distribution  at  Ti?  with  the  drift 
velocity  iteratively  updated  to  match  the  input 


T0=aTw  +  (l-a)Ti 


plenum  pressure  at  the  boundary.  Boundary  CD 
is  a  plane  of  symmetry  (assuming  a  multi-slot 
configuration).  The  farfield  vacuum  boundary 
DJ  is  placed  approximately  10  slot  widths 
downstream  of  the  slot  exit.  Each  wall  surface  is 
a  reflector  defined  by  an  accommodation 
coefficient  a  and  an  input  temperature,  either  Tj 
or  Tw  (surface  FI  only). 

The  key  feature  that  enhances  bulk  flow  in  the 
propellant  gas  is  the  differential  heating  of  the 
thin  film  heater  FI  (at  Tw)  relative  to  that  of  the 
nominal  plenum  inlet  and  surface  temperature  (at 
Tj).  In  this  study,  all  surfaces  except  the  thin 
film  heater  are  assumed  to  be  at  the  gas  inlet 
temperature  T^  More  accurate  simulations, 
which  iteratively  couple  the  gas  and  structural 
heat  transfer  properties  and  thus  eliminate  the 
need  for  an  input  surface  temperature 
distribution,  will  be  reported  at  a  later  time. 

The  code  is  instrumented  to  directly  sample  and 
spatially  resolve  the  flux  of  molecular  mass, 
momentum,  and  energy  through  an  arbitrary 
flowfield  plane,  here  typically  chosen  to  be  the 
slot  throat  (GG’)  or  the  slot  exit  (HIT).  The 
code  also  samples  the  axial  and  transverse 
velocity  distribution  function  at  arbitrary 
locations,  including  surface  boundaries,  to  allow 
more  detailed  comparison  with  the  simple  free 
molecule  theory.  In  this  manner,  performance 
losses  can  be  attributed  to  various  components  of 
the  geometry  such  as  the  slot  walls  GH. 

For  the  present  highly  rarefied  operating 
conditions,  the  DSMC  requirement  that  the 
nominal  cell  size  be  on  the  order  of  the  local 
molecular  mean  free  path  and  that  the  simulation 
timestep  be  small  relative  to  the  mean  collsion 
time  are  easily  met.  The  nominal  simulation  grid 
contained  2000  cells,  with  cell  sized  to  be 
typically  much  less  than  25%  of  the  local  mean 
free  path.  The  nominal  timestep  was 
approximately  50%  of  the  collision  frequency. 

The  simulations  were  run  for  2000  unsteady  and 
104  additional  steady-state  timesteps  during 
which  sampling  of  the  flowfield  occurred.  The 
nominal  simulation  contained  2  x  105  particles  at 
steady  state.  Statistical  scatter  in  the  predicted 
specific  impulse  is  estimated  to  be  much  less 
than  one  percent  due  to  the  large  sample  sizes 
obtained. 

For  simplicity  in  these  initial  design  studies, 
calculations  have  been  made  assuming  an  argon 


propellant,  although  it  is  far  from  an  optimum 
choice  for  an  electrothermal  thruster.  The 
typical  propellant  inlet  temperature  and  pressure 
are  300  K  and  50  Pa,  respectively.  The  nominal 
FMMR  geometry  has  a  slot  width  of  w  =  100 
pm,  an  expansion  angle  a  —  and  an  ^ 
expansion  slot  thickness  of  250  pm.  The  inlet 
conditions  and  slot  width  yield  a  Knudsen 
number  of  unity. 

GAS-SURFACE  INTERACTION 
MODELING 

Despite  advances  in  micromachining  technology, 
silicon  and  silicon-based  materials  are  still  the 
most  common  material  used  in  MEMS 
fabrication.  Exposure  of  pure  silicon  to  potential 
propellant  gases  such  as  ammonia  (NH3)  or 
water  (H20)  will  lead  to  dissociative  adsorption 
and  subsequent  nitridation  or  oxidation  of  the 
silicon  surface  even  under  benign  conditions 
(e.gftpbom  temperature).4  Very  little  information 
is  available  regarding  their  influence  on 
macroscopic  processes  such  as  gas  transport.5 
Molecular-level  data  indicate  desorption  is,  in 
general,  a  highly  non-equilibrium  process,6 
further  complicating  model  development. 

As  such  for  the  present  performance  studies,  an 
inert  surface  material  defined  by  a  wall 
temperature  Tw  or  T*  and  a  single 
accommodation  coefficient  a  are  assumed.  An 
estimate  of  the  influence  of  the  reflection  (or 
desorption)  process  on  overall  performance  is 
then  made  by  comparative  calculations  for  three 
models:  the  Maxwell  model,3  and  two  forms  of 
the  Cercignani-Lampis-Lord  (CLL)  model.7,8 

In  the  Maxwell  model,  a  fraction  a  of  incident 
particles  fully  accommodate  to  the  wall 
temperature  and  are  diffusely  reflected  with  the 
remainder  being  specularly  reflected.  Molecular 
beam  studies  have  shown  that  this  model  is  not 
realistic,  and  that,  in  general,  reflected  molecules 
at  a  surface  exhibit  a  lobular  distribution  in 
direction  and  a  continuous  spread  in  energy. 

In  the  original  form  of  the  CLL  model,  reflection 
varies  from  specular  (a  =  0)  to  fully 
accommodating  (a  =  1)  as  with  the  Maxwell 
model.7  However,  the  CLL  original  model 
allows  for  separate  normal  and  tangential  energy 
accommodation  coefficients  (which  for 
simplicity  are  equivalent  in  this  study).  This 


i 


3 


form  of  the  CLL  model  produces  physically 
reasonable  distributions  of  direction  and  energy 
for  reflected  molecules  although  it  can  not  be 
used  to  describe  discrete  internal  energy  states  of 
the  molecule. 

A  generalized  form  provides  for  diffuse 
molecular  reflection  at  a  surface  with  incomplete 
energy  accommodation.8  This  CLL  diffuse 
model  retains  the  energy  distribution  obtained 
from  the  original  CLL  model  but  chooses  new 
polar  and  azimuthal  scattering  angles 
corresponding  to  diffuse  reflection.  Both  the 
CLL  original  and  CLL  diffuse  models  are 
compared  to  the  simple  Maxwell  model  to  assess 
the  predicted  FMMR  performance. 

In  principle,  the  details  of  the  energy  deposition 
and  desorption  process  can  be  directly 
incorporated  into  the  DSMC  gas-surface 
interaction  procedure  as  they  become  available. 
For  the  present  range  of  FMMR  conditions,  gas- 
phase  processes  such  as  recombination  are 
negligible  due  to  the  low  operating  pressures. 
Therefore,  the  key  parameters  required  for  more 
detailed  studies  are  the  net  velocities  or 
translational  energy  components  of  the 
(potentially  dissociated)  products  of  the 
desorbing  propellant  molecule,  rather  than  the 
specifics  of  the  pathways  by  which  the 
adsorption  and  decomposition  occur. 


element  is  a  deposited  platinum  surface  with  a 
thickness  of  approximately  2000  A. 

The  temperature  dependent  accommodation 
coefficient  depends  on  several  factors  including 
the  physical  characteristics  of  the  gas,  the  solid 
surface,  the  surface  temperature,  and  the  gas 
temperature.  From  previous  research,  the 
accommodation  coefficient  for  argon  impacting  a 
platinum  surface  ranges  from  0.719  to  0.911  at  Tw 
=  300  K.  At  the  FMMR  operating  temperature 
of  Tw  =  600  K,  the  accommodation  coefficient 
ranges  from  0.510  to  0.63. 11  There  is  no  available 
data  for  the  accommodation  coefficients  of  argon 
colliding  with  a  silicon  surface  which  is  required 
to  estimate  the  energy  losses  on  the  expansion 
walls. 

Due  to  uncertainties  in  the  measurements  and  in 
some  cases  a  complete  lack  of  data,  the 
calculations  presented  in  the  following  section 
include  full  parameteric  analysis  with  the 
accommodation  coefficient.  For  simplicity,  the 
accommodation  coefficients  of  a  given  case  are 
assumed  to  be  equivalent  for  both  the  thin  film 
heater  and  the  expansion  walls.  In  light  of  these 
difficulties,  the  experimentally  derived 
accommodation  coefficients  can  be  used  as  a 
general  guideline  in  assessing  predicted  FMMR 
performance. 


THERMAL  ACCOMMODATION 
OF  GASES  ON  SURFACES 


An  extensive  literature  exists  detailing  the 
measurements  of  the  thermal  accommodation 
coefficient  made  for  various  gas-surface 
systems.9'12  Although  there  have  been  a  number 
of  notable  measurements  made  for  monatomic  or 
simple  diatomic  gases  interacting  with  pristine 
metal  surfaces,  obtaining  measurements  and 
developing  numerical  models  for  complex  flows 
on  engineering  surfaces  is  an  area  of  active 
research.13'16 


For  the  analysis  presented  here,  the  FMMR  is 
operating  on  an  argon  propellant  for  simplicity. 
With  current  MEMS  fabrication  practices,  the 
thin  film  heating  element  can  either  be  made 
from  a  polysilicon  material  or  sputter  deposited 
metals  such  as  platinum  or  gold.  It  is  assumed 
for  the  purposes  of  this  research  that  the  heating 


RESULTS 


Baseline  estimates  of  the  FMMR  performance 
(i.e.  thrust  and  specific  impulse)  and  flowfield 
behavior  were  obtained  assuming  a  Maxwell 
^  surface  model  with  complete  accommodation  (a 
^2  0T).  At  this  condition,  the  Maxwell  and  both 
%  vernons  of  the  CLL  model  are  equivalent. 
S^Fift&5-3  shows  the  FMMR  flowfield  contours 
for  complete  surface  accommodation  with  Tw  = 
600  K.  The  contours  consist  of  raw  data  with 
translational  temperature  on  the  left  side  and 
axial  velocity  on  the  right  side.  At  this  level  of 
rarefaction  (Kn  ~  1),  slip  phenomena  is  expected 
to  be  significant  as  confirmed  by  the  peak 
temperature  in  the  gas  near  the  heated  surface 
remaining  much  lower  than  Tw.  Figure  4 
compares  the  calculated  specific  impulse  with 
the  theoretical  value  (Eq.  4)  as  a  function  of  the 
heated  wall  temperature.  For  the  nominal  case, 
the  corresponding  thrust  per  unit  expansion  slot 
length  is  nearly  2.5  mN/m. 
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Figure  5  shows  the  centerline  axial  velocity 
distribution  function  at  the  expansion  slot  exit 
(H*  in  Fig.  2)  and  throat  (G’  in  Fig.  2)  for  the 
nominal  case  using  the  Maxwell  model  with  a  = 

1.  Also  shown  are  the  theoretical  half- 
Maxwellian  free  molecule  distributions 
evaluated  at  Tw  and  TV  The  theoretical  free 
molecule  axial  velocity  distribution  function 
would  consist  of  a  half-Maxwellian  at  Tw. 
Obviously,  the  flowfield  at  the  exit  differs  from 
the  simple  free  molegule  analysis;  however,  the 
comparisons  in  Fig£  4  suggest  that  they  can  serve 
well  for  initial  performance  studies. 

Figure  6  shows  the  net  axial  force  along  the 
expansion  slot  walls  for  various  expansion 
angles  a.  For  these  results,  the  heated  wall 
temperature  is  Tw  =  600  K  and  the  expansion  slot 
wall  temperature  is  =  300  K  with  full  surface 
accommodation.  The  net  increase  or  decrease  in 
thrust  contributed  by  the  expansion  walls  is 
given  by  the  integrals  under  the  curve  in  the 
figure.  For  smaller  expansion  angles  (a  =  40°  or 
54.74°)5  the  expansion  walls  contribute  a  net 
positive  increment  to  performance  by  limiting 
divergence  losses.  It  is  interesting  to  note  that 
the  enhancement  reaches  a  maximum  midway 
along  the  expansion  indicating  that  thin 
expansion  slots  are  preferred.  As  the  expansion 
angle  increases  for  a  fixed  slot  thickness,  the 
shear  forces  begin  to  dominate  and  can 
significantly  degrade  the  thruster's  performance. 
From  this  analysis,  relatively  thin  expansion 
walls  with  expansion  angles  a  <  60°  are 

desirable. 

As  mentioned  earlier,  the  characteristics  of  the 
FMMR  flowfield  and  performance  are  expected 
to  depend  on  the  details  of  the  gas-surface 
interaction  process.  Figure  7  shows  the  specific 
impulse  of  the  FMMR  as  a  function  of  the 
surface  accommodation  coefficient  for  the 
Maxwell  and  the  two  versions  of  the  CLL  model. 
The  models  differ  by  only  a  few  percent  at  a 
given  accommodation.  The  accommodation 
coefficient  in  this  case  is  assumed  to  be 
independent  of  the  heated  surface  temperature; 
with  this  assumption,  the  peak  specific  impulse 
occurs  for  a  =  1.  The  difference  between  the 
results  from  the  DSMC  incorporated  models  and 
the  free  molecule  theory  lies  in  the  effects  of  the 
expansion  slot  walls  on  the  FMMR  performance. 
As  indicated  in  Fig.  6,  the  expansion  slot  walls 
act  to  improve  the  performance  for  small 
expansion  angles.  Since  the  free  molecule 


theory  is  only  truly  valid  for  an  infinitely  thin 
orifice  expansion,  the  effect  of  the  slot  walls  is 
not  accounted  for  in  the  theoretical  results  shown 

Normalized  specific  impulse  profiles  across  the 
slot  expansion  exit  plane  (HH’  in  Fig.  2)  for 
various  conditions  are  shown  in  Fip T.  As 
expected,  the  maximum  specific  impulse 
contribution  is  on  the  FMMR  centerline  where 
losses  due  to  the  transverse  velocity  and 
collisions  with  the  expansion  walls  are 
minimized.  The  specific  impulse  drops 
approximately  30%  from  the  centerline  value  as 
the  expansion  walls  are  approached. 
Interestingly,  it  appears  that  the  nature  of  the 
scattering  process  has  a  larger  affect  on  the 
profile  shape  that  the  magnitude  of  the 
accommodation. 

Figure  9  shows  the  influence  of  the  model  on  the 
net  axial  force  along  the  expansion  slot  wall  for 
a  =  40°.  This  effect  is  larger  than  that  due  to 
modest  variations  in  the  expansion  sl£t  angle 
about  the  nominal  value  shown  in  Fi^  6.  Figure 
9  also  indicates  that  the  scattering  process 
contributes  significantly  to  the  determination  of 
the  optimal  slot  expansion  length.  To  minimize 
shear  and  heat  transfer  to  the  expansion  slot,  thin 
expansion  slots  are  desired.  However,  the 
benefit  of  re-directing  errant  molecules  in  the 
useful  thrust  direction  is  improved  for  longer 
slots. 

The  FMMR  flowfield  temperature  contours  for 
full  accommodation  and  for  the  CLL  £iffuse 
model  with  a  =  0.5  are  shown  in  FigJ/i 10.  As 
expected,  the  temperature  profile  near  the  heated 
wall  surface  shows  less  energy  is  transferred  into 
the  propellant  gas  under  incomplete 
accommodation  conditions. 

The  major  performance  loss  mechanism  is 
expected  to  be  due  to  heated  propellant 
molecules  colliding  with  relatively  cold  thruster 
walls,  particularly  in  the  expansion  region. 
Figure  11  shows  the  propellant  heat  flux  as  a 
function  of  axial  distance  along  the  FMMR 
expansion  slot.  The  entrance  of  the  slot  (i .e<jA 
position  closest  to  the  heated  wall)  is  denoted  as^ 
x/w  =  0.  Figure  11  shows  that  there  is  little 
deviation  between  the  gas-surface  interaction 
models  for  a  given  accommodation  indicating 
that  the  heat  flux  is  relatively  insensitive  to  the 
scattering  model  assumption.  As  expected,  the 
heat  flux  to  the  expansion  walls  varies  linearly 
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with  accommodation  coefficient  as  shown  in  Fig/ 
12.  For  example,  the  numerical  solution  for  a  = 
0.5  shows  approximately  one-half  of  the  heat 
flux  that  a  =  1.0  indicates. 


DISCUSSION 

As  shown  in  Figj  7,  the  details  of  the  gas-surface 
interaction  are  important  in  effectively 
determining  the  predicted  performance  of  the 
FMMR.  As  more  appropriate  propellants  are 
investigated  such  as  water  vapor,1  the  physics  of 
surface  interactions  will  become  increasingly 
complicated  by  adsorption,  dissociation  and 
internal  energy  state  excitation.  Currently, 
adequate  models  are  not  available  to  address 
these  interactions. 


The  mass  of  propellant  required  to  perform  a 
series  of  orbital  maneuvers  is  given  by 


mp=mi 


(  Av  'N 

IjpSo  _  | 

V 


(5) 


where  n^  is  the  dry  spacecraft  mass  (i.e.  without 
propellant)  and  Av  is  the  change  in  velocity  that 
the  propulsive  maneuvers  afford.  A  Av  of 
approximately  100  m/sec  is  typical  for  a 
nanospacecraft  mission  (n^  =  10  kg)  where  the 
micropropulsion  system  is  used  for  attitude 
control  over  the  mission  lifetime.  Figure  13 
shows  the  propellant  mass  required  for  the 
assumption  above  over  the  range  of  FMMR 
predicted  Isp  (Fig.  7). 

As  Eq.  -^/indicates,  an  error  in  estimating  the 
specific  impulse  can  lead  to  larger  errors  in 
predicting  the  required  propellant  mass  for  a 
spacecraft  mission  due  to  the  exponential  term. 
For  this  reason,  improved  gas-surface  interaction 
models  and  validating  data  are  required  to 
accurately  estimate  mission  parameters. 


assumed  than  the  details  of  the  surface 
scattering.  The  total  Isp  varies  by  approximately 
20%  over  range  of  accommodation  coefficients 
from  specular  (a  =  0)  to  diffuse  (a  =  1)  surface 
scattering.  However,  there  was  only  a  maximum 
difference  of  about  5%  between  the  models  for  a 
given  accommodation  coefficient.  Other  more 
microscopic  parameters,  such  as  axial  velocity 
distribution  functions,  appear  to  depend  more  on 
the  scattering  model  assumption. 

The  performance  of  the  FMMR,  for  a  given 
accommodation  on  all  surfaces,  is  maximized  for 
a  =  1  as  shown  by  the  numerical  results. 
Performance  may  be  further  optimized  by  using 
a  thin  film  heater  material  that  fully 
accommodates  while  choosing  a  more  specular 
surface  for  the  expansion  slots.  This  will  allow 
for  maximum  energy  input  into  the  gas  from  the 
heating  element  and  minimize  the  energy  lost  on 
the  expansion  walls. 

Additional  data  is  needed  to  assess  the 
accommodation  coefficient  of  gases  on  MEMS 
engineering  surfaces  such  as  silicon  (single  and 
polycrystalline),  silicon  oxide,  silicon  nitride, 
silicon  carbide,  and  certain  polymers.  More 
physical  gas-surface  interaction  models  are 
required  for  FMMR  calculations  which  simulate 
ideal  propellants  (typically  stored  as  liquid) 
interacting  with  surfaces  at  elevated  temperature. 
The  details  of  molecular  scattering,  adsorption 
(physi-  and  chemi-),  desorption,  dissociation  and 
internal  energy  mode  coupling  must  be 
addressed. 
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